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ABSTRACT 

X-ray grating observations have revealed great detail in the spectra of Novae in the Super Soft 
Source (SSS) phase. Notable features in the SSS spectra are blue-shifted absorption lines, P-Cygni 
line profiles, and the absence of strong ionization edges, all of which are indicators of an expanding 
atmosphere. 

We present, and make publicly available, a s et of 672 wind-type synthet ic spectra, obtained from 
the expanding NLTE SSS models introduced in I Van Rossum fc Nessl (j2010f ) with the PHOENIX stellar 
atmosphere code. The set presented in this paper is limited to solar abundances with the aim to 
focus on the basic model parameters and their effect on the spectra, providing the basis upon which 
abundance effects can be studied using a much bigger non-solar set in the next paper in this series. 

We fit the wind-type spectra to the five grating spectra taken in the SSS phase of nova V4743 
Sgr 2003 as an example application of the wind-type models. Within the limits of solar abundances 
we demonstrate that the following parameters are constrained by the data (in order of decreasing 
accuracy): column density Nn, bolometric luminosity Lboi, effective temperature T c g, white dwarf 
radius R, wind asymptotic velocity Dqo, and the mass- loss rate M. The models are also sensitive to 
the assumed white dwarf mass Mwd but the effect on the spectra can largely be compensated by 
the other model parameters. The wind-type spectra with solar abundances fit the data better than 
abundance optimized hydrostatic models. 

Subject headings: radiative transfer - novae, cataclysmic variables - stars: individual (V4 743 Sgr) - 
X-rays: stars 



1. INTRODUCTION 

Novae are cataclysmic variables caused by a white 
dwarfs in a binary systems that experience a thermonu- 
clear runaway in a compact surface-layer of hydrogen 
that was accreted from the companion. The outburst is 
detected as a strong sudden increase in optical brightness 
followed by a gradual decline over a time scale of months. 
Nova light curves vary s ubstantially from system to sys- 
tem (|Strope et al.l 12010). Some novae go through a Su- 
per Soft Source (SSS) phase, in which spectra are mea- 
sured similar to those of Super Sof t X-ray binary Sources 
ijKahabka fc van den Heuvel|[l997ft . 

Since the advent of the X-ray grating observato- 
ries, Chandra and XMM-Newton, 'high'-resolution X- 
ray spectra have been obtained for a number of novae 
in their SSS phase. These spectra contain a wealth of 
detail which have provided evidence of continued mass 
loss during the SSS phase: blue-shifted absorption lines, 
the abse nce of sharp ioni zation edges, and P-Cygni line 
profiles (|Ness et alj|2007[ ). 

These data so far were interpre ted with the help 
of hydrostatic atmo s phere models dRauch et al.l 120101 : 
lOsborne et all l20Tlt INess et al.l MlW . These models 
have proven to be a powerful tool for determining effec- 
tive temperatures and chemical compositions for grat- 
ing observations. But in order to match observed 
blue-shifted absorption lines the sy nthetic spectra are 
blue-shifted in an ad-hoc fashion (jNelson et al.l 120081 : 
iRauch et alll2010t INess et al.1 feOlll h which is not phys- 
ically consistent with having a velocity field in a radi- 



ally extended (vs. a compact hydrostatic) atmosphere. 
The effects of the velocity field on the physical condi- 
tions of the atmosphere (e.g. density and temperature 
stratification) can significantly change the spectrum and 
affect the interpretat i on of the observations as shown by 
IVan Rossum fc Nessl (pOlOf) . 

IPetz et al"n f2005) introduced expanding atmosphere 
model s for the SSS phase calculated with the PHOENIX 
codejHauschildt et al.l 119921: lHauschildt fc Baronlll99l 
I2004fh Those models are similar to the PHOENIX mod- 
els used to simulate ea rly UV and optical nova spectra 
(|Hauschildt et al.lll992D but applied to the very high tem- 
perature conditions of SSS. Although these models natu- 
rally reproduce blue-shifted absorption lines, neither the 
lines nor the continuum of the synthetic spectra match 
the shape of the observations we ll. 

In IVan Rossum fc Nessl (|2010f) a wind-type (WT) ex- 
panding atmosphere model was introduced for X-ray no- 
vae using the same PHOENIX code. Those models were 
shown to provide reasonable fits to two different SSS no- 
vae (RS Oph and V2491 Cyg), but their models were 
called preliminary and have since then neither been ap- 
plied and examined, nor been available for the public to 
use for that purpose. 

In this paper a large set of WT spectra is presented 
(and made publicly available), their properties explored, 
and their use demonstrated by application to the five 
V4743 Sgr 2003 grating observations. V4743 Sgr 2003 
was thoroughly analyzed with hydrostatic atmospheres 

1 The abundances were tuned to individual observations. 
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in iRauch et al.l (|2010f) , so that the results presented here 
can be directly compared to previous work. 

2. THE "WIND-TYPE" (WT) SSS NOVA MODEL 
ATMOSPHERE SETUP 

The WT model atmospheres are based on two obser- 
vations. Firstly, (at least) a part of the atmosphere 
is expanding. Secondly, given the high temperatures 
and small eff ective radii typ i cally determined for SSS 
spectra (e.g. lOsborne et all (|2011h . iNess et al.l (|2011f ). 
I Van Rossum fc Nessl (|20ldl V). (at least) a part of the at- 
mosphere is very close to the compact surface of the white 
dwarf. Based on these two observations the WT model 
atmospheres assume a spherically symmetric compound 
structure with a hydrostatic base and an expanding en- 
velope. The hydrostatic base is described with the clas- 
sical parameters effective temperature T e gf, gravitational 
acceleration log(g), and radius R, using the hydrostatic 
equilibrium condition 

^f = -pg for r<R, (1) 
ar 

where r is the radial distance to the center of mass, R the 
outer radius of the static core, p is the gas pressure, and 
p the density of the gas. The envelope is parameterized 
by a constant mass-loss rate M and a beta-law velocity 
field (see e.g. lLamersI (|1997| )) 

v(r>R)=v + (v oo -v )(l--\ , (2) 

with Voo the asymptotic velocity, and Vq the velocity at 
the base of the wind. The parameter (3 determines the 
length-scale in which the terminal velocity is reached. 
The density of the envelope follows from the continuity 
equation for a constant mass loss rate 

M 

p(r >R)= . (3) 

47rr z v(r) 

The base velocity vq is refined in the model to provide 
a smooth run of the density and its gradient over t he 
transition boundary at r — R (see lVan Rossuml |2009)). 
This is not possible for any arbitrary value of /3. Ex- 
periment shows that f3 = 1.2 ± 0.05 for the WT models 
to have a smooth transition between static core and ex- 
panding envelopfl This is slightly higher than the value 
of P — 0. 9 ± 0.1 for radiat i on dr iven winds of blue giants 
see e.g. iPauldrach et~aH (|l~986r ib 

In conclusion, the model atmosphere is defined by five 
parameters^, the three classical hydrostatic parameters: 
T g, log(<?), R, and two parameters for the expanding 
envelope: M, v^. 

3. CHARACTERISTICS OF THE WT MODEL 

With five parameters, the WT models have much more 
freedom than static atmospheres (which have three) . The 

2 Note that the model atmosphere, and thus it's spectrum, 
depends on the assumed velocity law. The effect of (3 on the spectra 
is small within the allowed range. The effect of assuming a different 
type of velocity-law is potentially bigger. This is not explored 
within the scope of this work. 

4 The chemical composition adds a large number of additional 
free parameters to the model, the discussion of which is deferred 
to next paper of this series. 



Table 1 

Overview on the effects of the five model parameters on four 
observables. 
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Notes: This table serves as a guide in the process of finding good 
spectral fits to data, ~f means a positive correlation, J, a negative 
correlation, and — no correlation between the model parameter and 
the observable (abstract property of the spectrum). The double ar- 
rows ff denote the parameters that are most characteristic for each 
observable. 

The parameter log(g) has no double arrows meaning that it is not 
well constrained by the observables. The Color is affected by three 
parameters, which means that T ef j generally is not determined 
uniquely, but depends on the value for log(cj). 

effect of each of these p arameters on the m odel atmo- 
sphere was examined in iVan Rossuml ((2009) . Here we 
demonstrate the effect on the observables, i.e. the spec- 
tra, focusing on four abstract properties: 

1. Luminosity (or brightness) 

2. Color (or blueness) 

3. Blue-shift of absorption features 

4. Emission wing prominence in P Cyg profiles 

In the following, the details of the relations between 
model parameters and these observables is given, and 
an overview is shown in the Table [T] This table serves as 
a guide in the process of fitting WT spectra to data. 

Radius — The atmosphere is self-similar with variation 
of R, meaning that it does not affect the flux produced 
by the model. Yet, the total luminosity 

iboi = 4^i?V sb r c 4 ff , (4) 

with cr s b the Stefan-Boltzmann constant, scales trivially 
with R 2 . 

Effective Temperature — T e g makes the spectrum both 
brighter and bluer (or dimmer and redder) in a way sim- 
ilar to how a blackbody spectrum changes with increas- 
ing (or decreasing) temperature. But, in addition, the 
ionization balance of all chemical species changes with 
temperature which has a strong influence on the spectra. 
The reason is that with effective temperatures typical 
for SSS models the C, N and O atoms are almost com- 
pletely ionized to bare nuclei, which do not contribute 
to the opacity. If the ionization balance changes slightly 
from bare nuclei towards H-like C, N and O then bound- 
free transitions start to significantly contribute to blue 
opacities making the spectrum redder. 

The brightness variation with effective temperature 
can completely be compensated by the radius R, so that 
T e ff essentially is a color parameter. The color variation 
of spectra with T e g is shown in Figure [l] where spectra 
are scaled with T~ s to identical luminosities. 

Surface Gravity — Lowering \og{g) makes the atmosphere 
less compact, so the density gradient becomes smaller. 
As mentioned above, the wind base-velocity vq is refined 
in the model to provide a smooth run of the density and 
its gradient over the core-envelope boundary. A smaller 
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The color of the spectrum changes bluer with increased T c ff much 
more quickly than a blackbody does. This is caused by the high 
temperatures ionizing key elements like carbon, nitrogen and oxy- 
gen to bare nuclei, removing their bound-free opacities from the 
blue side of the spectrum if the temperature is increased even fur- 
ther. 

density gradient in the base of the envelope is obtained 
with higher v values (from Equations ((2]) and (jSJ), so 
that the transition density becomes lower. Consequently, 
a lower surface gravity gets a wind that is optically thin- 
ner, the hot compact core moves "up" on the optical 
depth scale, and the spectrum becomes bluer. The re- 
verse holds for increasing log (g). 

The assumed white dwarf mass follows from the values 
for log(g) and R as 

Mwd = ^ (5) 

where G is the gravitation constant. In Figure [2] the 
color variation of spectra with \og(g) is shown. Given the 
large spread in white dwarf masses the influence is rather 
small. Furthermore, as shown in Figure G2 the color effect 
of the surface gravity can largely be compensated by the 
color effect of the effective temperature, leaving log(p) 
poorly constrained by the observables. 

Asymptotic Velocity — linearly scales the velocity law, 
apart from a small offset vq <C . For a fixed mass- loss 
rate M the density of the wind changes with Voq. But 
the density profile of the wind has a much bigger influ- 
ence on the spectrum than the velocity profile. In order 
to compare different velocities with identical density pro- 
files the mass-loss rate is to be changed along with the 
terminal velocity, keeping M/voo fixed. In Figure [4] the 
effect on the blue-shift of absorption profiles is shown. 

Mass-loss Rate — M is proportional to the wind den- 
sity. The emission wings of P Cyg profiles originate 
from the radially extended atmosphere, from those parts 
that are not in the line of sight to the compact core. 
With increasing mass-loss rate the extended atmosphere 
becomes denser, therefore more material is available to 
line-scatter photons toward the observer, and so emission 
wings of P Cyg profiles become stronger. 
With more prominent emission wings, the absorption 




10 15 20 25 30 35 40 

Wavelength [A] 



Figure 2. The effect of log(g) on the WT spectra. A lower sur- 
face gravity in the model gets a matching wind that is optically 
thinner, as explained in the text. Consequently, the compact, hot, 
static part of the atmosphere moves "up" on the optical depth scale 
leading to a bluer color of the spectrum. 

Although this is a significant effect per se, the associated range 
in underlying white dwarf masses is very large (a factor of 2.2). 
Furthermore, the color effect of the surface gravity can largely be 
compensated by the color effect of the effective temperature as 
shown in Figure [3] This leaves log(g) poorly constrained by the 
observables. 
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Figure 3. The color effect of the surface gravity on the WT spec- 
tra (see Figure [2J, can largely be compensated by the color effect 
of the effective temperature (see Figure [TJ . The fluxes of each 
pair of synthetic spectra are scaled to identical luminosities with 
(700kK/T cff ) 4 (black-red) and (500kK/T eff ) 4 (green-blue) respec- 
tively, which corresponds to changing R. 

The underlying white-dwarf mass differs by a factor of 5.1 (black- 
red) and 2.7 (green-blue) respectively, log(g), and therefore M\yrj, 
are poorly constrained by the observables. 

wing of the line profiles may appear to blue-shift more 
strongly with increasing M. 

Another result of the increased wind density is that the 
transition to the compact core happens at higher den- 
sities and higher temperatures. The temperature and 
density gradients in the wind are lower than in the static 
core, so that with M these increase throughout the whole 
wind region. Both a higher temperature and a higher 
electron density boost ionization, shifting the balance 
towards bare nuclei and removing their opacity. Conse- 
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Wavelength [A] 

Figure 4. The effect of Voo on the WT spectra. In order to leave 
the density profile unaffected in the comparison M is scaled linearly 
with described in the text. 

The line profiles become broader and the absorption part more 
blue-shifted with increasing Uoo. Also, narrow features get more 
smeared out so that high-velocity spectra are smoother and low- 
velocity spectra show more small line features. 

quently, the spectra turn bluer in color as the mass-loss 
rate increases. 

All three effects M has on the spectra are shown in fig- 
ure The spectral color and the blue-shift can be partly 
compensated by T e g and u_o respectively. Therefore, the 
most characteristic effect of M is the prominence of emis- 
sion wings in P Cyg profiles. 

4. COMPARISON WITH OBSERVATIONS 
4.1. Interstellar absorption: NOT a free parameter 

The colors of SSS spectral observations are strongly 
affected by absorption from the interstellar material 
(ISM). In the comparison with observed data the syn- 
thetic flux is to be multiplied with the wavelength de- 
pendent transmission coefficient. Generally, this co- 
efficient depends on the ISM column depth Ah and 
the chemical composition of the ISM. In this paper 
we use the ISM absorption cross- s ection routines from 
IBalucinska-Church fc McCammonl ([1992) with solar (or 
almost-solaiQ) abundances. 

The ISM-absorbed spectrum is very sensitive to the Ah 
parameter, especially the slope of the red tail (the part 
redder than about the peak wavelength) of the spectrum. 
Although the red tail also is affected by the WT model 
parameters T e g, log(g), and M, all of those parameters 
have a bigger impact on the blue half of the spectrum. 
We find, as a consequence, that the value of Ah can 
accurately be determined, separately from the WT model 
parameters. However, this can only be done if the red tail 
is completely measured by the instrument, like with the 
LETGS spectrograph on board CHANDRA. The RGS 
spectrograph on board XMM- Newton is limited to 38. 2 A 
and therefore does not capture the red tail of the SSS 
spectrum, leaving Ah undetermined. 

5 For the V4743 Sgr 2003 data, carbon and nitrogen abun- 
dances of [C] = -1 and [N] = 0.4 are assumed in order to match the 
non-absorbed fluxes on the right-hand and left-hand side of either 
ionization edge. All other abundances are assumed to be solar. 



Fit procedures that treat Ah as a free parameter are 
sensitive to t he w avelength range, as shown e.g. in 
iRauch et all (|2010l) . Figure 12. That same Figure 12 
also shows that chi-squared methods do not accurately 
determine Nn even when the red tail is available, as 
evident from the non-matching slopes for wavelengths 
greater than 30A. The problem is that the red tail has 
very little counts so that chi-squared methods, which op- 
timize the fit in the high-counts bins, are not sensitive to 
an accurate slope in the red tail of the spectrum. For this 
reason we claim that Ah should not be treated like a free 
parameter. Doing so gives the fit too much freedom and 
leads to incorrect fit parameters. Instead, once Ah is de- 
termined from the long wavelength tail (provided this is 
possible for the data) , fitting the short wavelength range 
with Ah fixed will give more realistic fit parameters. 

The sensitivity of the ISM-absorbed spectrum to Ah is 
demonstrated in Figure |6l where a synthetic spectrum is 
compared with an observation using three different val- 
ues. In this example Ah can be determined with an 
uncertainty of ±3%. The associated uncertainties in R 
and Lboi are ±1.5% and ±2.5% respectively. 

4.2. Counts spectra, Flux spectra, or both 

Comparisons between observations and synthetic spec- 
tra are often presented in the "as-observed" representa- 
tion, as ISM-absorbed counts per bin on a linear scale, 
or in the "as-modeled" representation, as ISM-absorption 
corrected flux on a logarithmic scale. Both representa- 
tions have their advantages. 

In the "as-observed" representation synthetic spectra 
are 'forwardly' corrected for ISM-absorption and con- 
volved with the instrumental response of the detector. 
The counts per bin numbers indicate how well each bin 
was exposed, preserving the statistical nature of X-ray 
measurements. Note that any chi-squared fits in this 
representation rely on the arbitrariness of the detector 
properties and the ISM absorption. 

In the "as-modeled" representation observations are 
'reversely' corrected for ISM-absorption and the instru- 
mental response of the detector. Deconvolution of the in- 
strumental response is complicated and is compromised 
by small-number statistics. Given that the observation is 
exposed well enough this representation shows the true 
flux. Sure enough, a good model needs to not only repro- 
duce the arbitrary (instrumental response and ISM ab- 
sorption) high-count bins, but the whole run of the true 
flux from the source. Also, the logarithmic scale bet- 
ter highlights the dim blue end of the spectrum which is 
sensitive to details of the atmosphere (e.g. temperature, 
mass-loss and composition). 

For these reasons, we advocate to show comparisons in 
both the "as-observed" and the "as-modeled" represen- 
tations. 

5. THE WT MODEL PARAMETER GRID 

With this paper a moderately sized parameter grid of 
WT model spectra is introduced and made publicly avail- 
ably. It is a solid four-dimensional parameter grid cov- 
ering the four WT model parameters other than R (since 
variation with R is trivial as explained in Section [3]). For 

6 http://flash.uchicago.edu/~daan/WT_SSS/ 
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Figure 5. The effect of M on the WT spectra. As described in the text, three of the abstract observables are affected by the mass- loss 
rate. Firstly, the emission wings of P Cyg profiles get stronger. Secondly, the blue-shift of absorption features increases slightly. Thirdly, 
the color of spectra turns bluer. 

Since the spectral color and the blue-shift can be partly compensated by T e g and respectively, the most characteristic effect of M is 
the prominence of emission wings in P Cyg profiles. 

The fluxes are scaled to identical luminosities with a factor of (700kK/T e fj ) 4 . 




Wavelength [A] 

Figure 6. The ISM absorption parameter 7Vh strongly affects the comparison of synthetic spectra with data. The plot shows a V4743 
Sgr (day 302) count spectrum compared with three WT spectra using different values of the interstellar absorption Nn: 1-44 (red), 1.36 
(green), and 1.28 (blue) in units of 10 21 cm — 2 . The values of T cf j , v^, M, L ho i, and R shown in the plot are in the units: K, km/s, Mq/yi, 
erg/s, and km. 

Nft determines the slope of the red tail of the spectrum. In this case the observed red tail slope is reproduced with Nn = 1.36 ± 0.04 X 
10 21 cm -2 (green curve). 
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Table 2 

WT spectra are calculated on these parameter grid points. 



Par am. 


# 


Grid point values 


T c ff 


8 


[450, 475, 500, 550, 600, 650, 700, 750] xlO 3 K 


l°g(9eff) 


7 


-15.55 + T 4 ff + [0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2] 


Voo 


4 


[1200, 2400, 3600, 4800] km/s 


M 


3 


«oo/2400 X [HT 8 , 3 ■ 10~ 8 , 10~ 7 ] M /yr 



improved stability of the WT model log(g) is not directly 
specified but instead of g the effective gravitational ac- 
celeration geff = g — arad (at the surface of the static 
core), where a rac j is the acceleration from radiation pres- 
sure. Therefore, the parameter grid is spaced in terms of 
log(g c ff). Because the radiation pressure roughly scales 
with T^ s we also scale g e g with T e 4 ff . This way a gravita- 
tionally loosely bound atmosphere remains loosely bound 
along the T c g grid axis, i.e. the ratio g/a ra d is roughly 
independent of T e s. 

The parameter grid size is 8x7x4x3 in the dimensions 
T c g, log((7eff), Voo, and M respectively, for a total of 672 
spectra. The parameter values at the grid points are 
listed in Tabled 

6. UNCERTAINTIES FROM M WD -AMBIGUITY OF WT 
MODELS 

In Section [3] it was explained that log(g) affects the 
WT spectra (Figure [2]), but that the effect can be largely 
compensated with T e g (Figure [3]). But with T e g the lu- 
minosity Lboi of the model changes which can be com- 
pensated by changing R. Thus, the white dwarf mass 
changes with log(g) and R simultaneously (see Equation 
([5]) ) . The ambiguity of the white dwarf mass in the WT 
models is demonstrated in Figure[71 where a 0.47M Q and 
a 1.42M model give an equally good description of the 
data. The differences between the WT spectra are too 
small to favor one above the other, especially considering 
the enormous difference in white dwarf mass. 

The MwD-ambiguity imposes uncertainties in the cor- 
related quantities T e g and i?, the size of which is to be 
determined from fits to the data. From Figure[7]we deter- 
mine that the uncertainty from A/vvD-ambiguity is ±5% 
in T er f and about ±8% in R. The uncertainty from A/wd- 
ambiguity in the bolometric luminosity Lboi is approxi- 
mately ±1.5%. 

The assumption of solar abundances causes additional 
systematic uncertainties, the sizes of which have to be 
determined from fits of non-solar models to the data, 
which is beyond the scope of this work. 

7. FITS TO V4743 SGR 

Figure [8] and [9] demonstrate the determination of Vqc 
and M from the Day 302 spectrum of V4743 Sgr. 

With increasing wind velocity some absorption features 
get increasingly blue-shifted while other lines only change 
in profile shape, depending on where in the atmosphere 
the lines are formed. Lines that become optically thick 
in the outer wind layers are more blue-shifted than those 
formed deeper inside. In the high-velocity models the 
blue-shift in the lines that are affected is stronger than 
what is observed for V4743. Also, in those models the 
line profiles are too washed-out (or smeared-out), which 
is most apparent in the lines that do not blue-shift. 



Table 3 

Best fit parameters and uncertainties for all five grating 
observations of V4743 Sgr. 



Day 


N H 






R 


Voo 


M 


180 


1.34 


6.65 


5.5 


10. 


1.2 


1.6 


196 




5.88 


5.5 


9.5 


1.2 


1.6 


302 


1.36 


6.13 


5.5 


9.7 


1.2 


1.6 


371 


1.39 


3.77 


5.5 


7.6 


1.2 


1.6 


526 


1.31 


1.93 


4.8 


7.3 


1.2 


1.6 


uncert. 


±3% 


±3% 


±8% 


±11% 


±30% 


±50% 



Notes: 

The uncertainties do not include the effects from abundance 
variations, which are not explored in the scope of this work, nor 
the distance uncertainty. The assumed distance to V4743 Sgr is 
4kpc (see text). 

Nfi is in units of 10 21 cm -2 . ISM abundances are assumed to be 

solar, except [C]=-1.0 and [O]=0.4 (see text). 

L bo i is in units of 10 37 erg/s 

T e g is in units of 10 5 K 

R is in units of 10 8 cm 

Voo is in units of 10 8 cm/s 

M is in units of lO _8 M /yr 

If the mass-loss rate increases each absorption line gets 
a more prominent emission wing. The combined effect of 
many such P Cygni profiles leads to an increasing number 
of peaks and pits in the spectra, and steep slopes between 
them. At the same time, the spectra become bluer (see 
section|3j , which is compensated by lowering the effective 
temperature. In turn the white dwarf radius is increased 
to maintain the overall flux level, and the surface gravity 
changes in order to compare models with approximately 
equal white dwarf masses. The relative uncertainty in M 
is 50% and the associated additional uncertainties in T c ff 
and R are 5%. 

The fit results for all five grating spectra available for 
V4743 Sgr are shown in Figures HU1 to H2l The white 
dwarf masses (for which the WT models are degenerate, 
see section |6]) in each of the fits are chosen around an 
arbitrary value of Mwd = 1.1M© as close as possible 
within the available grid spacing (see section [SJ. 

The Day 526 observation is less well exposed than the 
first four observations. We re-binned the Day 526 data 
by a factor of two to suppress the noise. The Day 196 ob- 
servation was made with the RGS spectrograph on board 
XMM-Newton with an upper wavelength limit of 38. 2A. 
For this observation no individual determination of Ah 
is possible (as explained in section S]). Therefore, we as- 
sumed the average value from the other observations. 

The best-fit parameters and relative uncertainties are 
listed in Table [31 The dist ance to V4743 is not accu- 
rately known in literature. Niclboc k fc Schmidtobreickl 
(2003) determine an upper limit of 6kpc from optical de- 
cline rate considerations and a range of 0.4-1.2kpc using 
reference stars, but also point out that these values are 
uncertain. We assume a value of D = 4kpc as a "golden 
mean" guess. Note, that values for R and Lboi scale lin- 
early and quadratically with D. The uncertainties for 
R and Lboi as listed in Table |3J only reflect effects that 
can be directly determined with the models. For final 
values, the distance uncertainty has to be propagated 
accordingly. 

8. DISCUSSION AND CONCLUSIONS 
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Figure 7. The fit of two WT spectra with different log(g) vaiues shows that observed data can be equally well fit with very different 
white dwarf masses of 0.47Mq (green curve) and a 1.42-Mq (red curve). The comparison with data from V4743 Sgr (day 203) is shown in 
two representations, a logarithmic unabsorbed flux plot called "as-modeled" (upper-panel) and a linear ISM-absorbed counts plot called 
"as-observed" (lower panel). The values of T e jj, w^,, M, L bo i, ancl R shown in the plot are in the units: K, km/s, M^/yr, erg/s, and km. 
The MwD _amrj iguity causes uncertainties in the fit parameters T cf j and R, while the determined value of L bo i is little affected. 
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R=9700 
A/wd=1.07 
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V4743Sgr day 302, LETGS, AMH=1.36e21 
7eff=550k, log(a) = 8.18, Unf = 1 2, «dot=1.6e-08 
reffsSSOk, log(o) = 8.18, Hnf = 24, J/dot = 3. 1 e-08 
reffsSSOk, log(a) = 8.18, Knf = 36, J/dot = 4.7e-08 
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Figure 8. Three WT model spectra with different wind asymp- 
totic velocities compared to the V4743 Sgr Day 302 observation: 
Vao = 1200 km/s (red), 2400 km/s (green), and 3600 km/s (blue). 
Note that the mass-loss rate changes along with Voc, as explained 
in section \3\ and Figure [4] The units of the quantities shown are 
described in the caption of Figure 

Some absorption features become increasingly blue-shifted with 
Voo> e.g those around 24.2, 26.3, and 27.5 A, whereas other remain 
at their rest wavelengths, e.g. those around 25.8, 28.2, and 31.1 
A. Furthermore, the spectral features become increasingly washed- 
out, i.e. broader and shallower, with increasing Voo. 
The steeper spectral features and the moderate blue-shifts in the 
■u 00=1200 km/s model give the best match to the data. 

The representation of the data by the models is good, 
especially for the first four observations. Many, though 
not all, of the observed spectral features are (at least) 
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R= 10300 
A/wd=0.96 
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R=9700 
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Lbol = 4.17 
J?=8000 
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V4743Sgr day 302, LETGS, jVH = 1.36e21 
7"eff = 550k, log(g)=8.08, Hnf=12, Mdot = 5.0e-09 
reff = 550k, log(s) = 8.08, Vinf=12, Afdot= 1 ,6e-08 
reff = 550k, log(j) = 8.07, Hnf = 12, J/dot = 5.0e-08 
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Figure 9. Three WT model spectra with different mass-loss rates 
compared to the V4743 Sgr Day 302 observation: M = 5.0 ■ 
10~ 5 M /yr (red), 1.6 ■ 10"* M©/yr (green), and 5.0- 10" 9 M /yr 
(blue). The units of the quantities shown are described in the cap- 
tion of Figure [7] Note that the wind densities are relatively high 
given the low wind velocities, Vac = 1200 km/s. 
With increasing mass-loss rate each isolated absorption feature gets 
a more prominent emission wing so that in the interplay of many 
overlapping lines the spectrum shows deeper and narrower features. 
In comparison with the data, this effect is too strong for the high- 
est mass-loss rate model and too weak for the lowest. 
Within the available parameter grid spacing a value of M = 
1.6 • 10 — 8 M0/yr reproduces the observations best. The match 
in detail also depends on the chemical abundances which are not 
optimized in the scope of this work. 
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Wavelength [A] 



Figure 10. WT model spectra fit to V4743 Sgr grating data in the "as-observed" representation (see section|4]l. The values of T e g, iioo, 
M, ibob an d R shown in the plot are in the units: K, km/s, Mq/jt, erg/s, and km. 
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Figure 11. WT model spectra fit to V4743 Sgr grating data in the "as-modeled" representation (see sectional. The units of the quantities 
shown are described in the caption of Figure ITOl 
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7eff = 550k, log(g) = 8.18, finf=12, Mdot=1.6e-08 



tronsmission — — — — 




V4743Sgr day 371, LETGS, AM = 1.39e21 Z 

reff = 550k, log(g) = 8.44, Hnf=12, Mdot=1.6e-08 - 

transmission — — — — 




Wavelength [A] 

Figure 12. Like the previous figure, but with a smaller wavelength range, zooming in on the detailed line profiles in the spectra. 
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roughly reproduced. The overall flux level matches well 
over the whole wavelength range except shortward of the 
22.5 A where the N VI ionization edge is not strong 
enough in the models. Shortward of 18.6 A, the NVII 
ionization edge, the discrepancy grows even stronger. 
Also, the strong NVI and NVII lines at 28.9 and 24.8 A 
are too weak in the models. Obviously, a super-solar N 
abundance would improve these aspects of the fits to the 
data. Note that the relative strengths of ionization edges 
and absorption lines are largely determined by NLTE ef- 
fects. The dependence of these on abundance variations 
is generally hard to predict. 

The ISM column depth varies as little as 6% over the 
four LETGS observations. 

The first four observations are best-fit by a single ef- 
fective temperature. This would not have been the case 
if the log(g) parameter would not have been adapted to 
satisfy a fixed underlying white dwarf radius. The varia- 
tions in bolometric luminosity can here be explained by 
only varying the white dwarf radius and adapting the 
surface gravity correspondingly. 

The mass-loss rate and the asymptotic velocity of the 
wind do not significantly vary over the five grating ob- 
servations of V4743 Sgr 2003. 

Overall, the WT models with solar abundances rep- 
resent the V4743 Sgr data better than current hydro- 
static models with optimized abundances. Though abu- 
dance optimizations will certainly further improve the 
fits, strong departures as determined with hydrostatic 
models are possibly not necessary. 

With five parameters the WT models have more "free- 
dom" than hydrostatic models, which have three. The 
two extra parameters make the procedure for finding 
good fits slightly more complex, especially since their 
effects on the spectra partially overlap. But we have pre- 
sented a straight forward way to uniquely determine the 
best-fit model parameters, apart from the white-dwarf- 
mass-ambiguity in the WT models that leave -Mwd un- 
determined. 

Comparisons between models and SSS data should 
preferably be done in both "as-observed" and "as- 
modeled" representations, because each highlight differ- 
ent important aspects of the spectra. Fits that only look 
good in one or the other representation are not accept- 



able. 

The two properties most accurately constrained by the 
data are: 1) the ISM column depth, and 2) the bolo- 
metric luminosity (provided the distance to the source is 
known accurately). The column depth can only be ac- 
curately determined for LETGS data, because the RGS 
wavelength coverage is insufficient. We claim that precise 
determination of the ISM column depth is very impor- 
tant (or even crucial) if accurate fit parameters are to be 
obtained. The additional freedom introduced by treating 
iV H like a free parameter can easily double or triple the 
uncertainties in most model fit parameters. 



We thank Sumner Starrfield for the helpful discus- 
sions and comments on the draft of this manuscript, Jan- 
Uwe Ness for providing the observational data, and John 
Dombeck for providing his IDL color tables. 
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